Deciphering "Immaturity-Stemness" in Human Epidermal Stem Cells at the Levels of Protein-Coding and Non-Coding Genomes: A Prospective Computational Approach.

The epidermis hosts populations of epithelial stem cells endowed with well-documented renewal and regenerative functions. This tissue thus constitutes a model for exploring the molecular characteristics of stem cells, which remain to date partially characterized at the molecular level in human skin. Our group has investigated the regulatory functions of the KLF4/TGFB1 and the MAD4/MAX/MYC signaling pathways in the control of the immaturity-stemness versus differentiation fate of keratinocyte stem and precursor cells from human interfollicular epidermis. We described that down-modulation of either KLF4 or MXD4/MAD4 using RNA interference tools promoted an augmented stemness cellular status; an effect which was associated with significant transcriptional changes, as assessed by RNA-sequencing. Here, we have implemented a computational approach aimed at integrating the level of the coding genome, comprising the transcripts encoding conventional proteins, and the non-coding genome, with a focus on long non-coding RNAs (lncRNAs). In addition, datasets of micro-RNAs (miRNAs) with validated functions were interrogated in view of identifying miRNAs that could make the link between protein-coding and non-coding transcripts. Putative regulons comprising both coding and long non-coding transcripts were built, which are expected to contain original pro-stemness candidate effectors available for functional validation approaches. In summary, interpretation of our basic functional data together with in silico biomodeling gave rise to a prospective picture of the complex constellation of transcripts regulating the keratinocyte stemness status.

Global miRNA expression reveals novel nuclear and mitochondrial interactions in Type 1 diabetes mellitus.

Background: Considering the potential role of miRNAs as biomarkers and their interaction with both nuclear and mitochondrial genes, we investigated the miRNA expression profile in type 1 diabetes (T1DM) patients, including the pathways in which they are involved considering both nuclear and mitochondrial functions. Methods: We analyzed samples of T1DM patients and control individuals (normal glucose tolerance) by high throughput miRNA sequencing (miRNome). Next, five miRNAs - hsa-miR-26b-5p, hsa-let-7i-5p, hsa-miR-143-3p, hsa-miR-501-3p and hsa-miR-100-5p - were validated by RT-qPCR. The identification of target genes was extracted from miRTarBase and mitoXplorer database. We also performed receiver operating characteristic (ROC) curves and miRNAs that had an AUC > 0.85 were considered potential biomarkers. Results: Overall, 41 miRNAs were differentially expressed in T1DM patients compared to control. Hsa-miR-21-5p had the highest number of predicted target genes and was associated with several pathways, including insulin signaling and apoptosis. 34.1% (14/41) of the differentially expressed miRNAs also targeted mitochondrial genes, and 80.5% (33/41) of them targeted nuclear genes involved in the mitochondrial metabolism. All five validated miRNAs were upregulated in T1DM. Among them, hsa-miR-26b-5p showed AUC>0.85, being suggested as potential biomarker to T1DM. Conclusion: Our results demonstrated 41 DE miRNAs that had a great accuracy in discriminating T1DM and control group. Furthermore, we demonstrate the influence of these miRNAs on numerous metabolic pathways, including mitochondrial metabolism. Hsa-miR-26b-5p and hsa-miR-21-5p were highlighted in our results, possibly acting on nuclear and mitochondrial dysfunction and, subsequently, T1DM dysregulation.

CircRNAs as Potential Blood Biomarkers and Key Elements in Regulatory Networks in Gastric Cancer.

Gastric cancer (GC) is the fifth most common type of cancer and the third leading cause of cancer death in the world. It is a disease that encompasses a variety of molecular alterations, including in non-coding RNAs such as circular RNAs (circRNAs). In the present study, we investigated hsa_circ_0000211, hsa_circ_0000284, hsa_circ_0000524, hsa_circ_0001136 and hsa_circ_0004771 expression profiles using RT-qPCR in 71 gastric tissue samples from GC patients (tumor and tumor-adjacent samples) and volunteers without cancer. In order to investigate the suitability of circRNAs as minimally invasive biomarkers, we also evaluated their expression profile through RT-qPCR in peripheral blood samples from patients with and without GC (n = 41). We also investigated the predicted interactions between circRNA-miRNA-mRNA and circRNA-RBP using the KEGG and Reactome databases. Overall, our results showed that hsa_circ_0000211, hsa_circ_0000284 and hsa_circ_0004771 presented equivalent expression profiles when analyzed by different methods (RNA-Seq and RT-qPCR) and different types of samples (tissue and blood). Further, functional enrichment results identified important signaling pathways related to GC. Thus, our data support the consideration of circRNAs as new, minimally invasive biomarkers capable of aiding in the diagnosis of GC and with great potential to be applied in clinical practice.

Can miRNA Indicate Risk of Illness after Continuous Exposure to M. tuberculosis?

The role of regulatory elements such as small ncRNAs and their mechanisms are poorly understood in infectious diseases. Tuberculosis is one of the oldest infectious diseases of humans and it is still a challenge to prevent and treat. Control of the infection, as well as its diagnosis, are still complex and current treatments used are linked to several side effects. This study aimed to identify possible biomarkers for tuberculosis by applying NGS techniques to obtain global miRNA expression profiles from 22 blood samples of infected patients with tuberculosis (n = 9), their respective healthy physicians (n = 6) and external healthy individuals as controls (n = 7). Samples were run through a pipeline consisting of differential expression, target genes, gene set enrichment and miRNA-gene network analyses. We observed 153 altered miRNAs, among which only three DEmiRNAs (hsa-let-7g-5p, hsa-miR-486-3p and hsa-miR-4732-5p) were found between the investigated patients and their respective physicians. These DEmiRNAs are suggested to play an important role in granuloma regulation and their immune physiopathology. Our results indicate that miRNAs may be involved in immune modulation by regulating gene expression in cells of the immune system. Our findings encourage the application of miRNAs as potential biomarkers for tuberculosis.

Comprehensive analysis of germline mutations in northern Brazil: a panel of 16 genes for hereditary cancer-predisposing syndrome investigation.

BACKGROUND: Next generation sequencing (NGS) has been a handy tool in clinical practice, mainly due to its efficiency and cost-effectiveness. It has been widely used in genetic diagnosis of several inherited diseases, and, in clinical oncology, it may enhance the discovery of new susceptibility genes and enable individualized care of cancer patients. In this context, we explored a pan-cancer panel in the investigation of germline variants in Brazilian patients presenting clinical criteria for hereditary cancer syndromes or familial history. METHODS: Seventy-one individuals diagnosed or with familial history of hereditary cancer syndromes were submitted to custom pan-cancer panel including 16 high and moderate penetrance genes previously associated with hereditary cancer syndromes (APC, BRCA1, BRCA2, CDH1, CDKN2A, CHEK2, MSH2, MSH6, MUTYH, PTEN, RB1, RET, TP53, VHL, XPA and XPC). All pathogenic variants were validated by Sanger sequencing. RESULTS: We identified a total of eight pathogenic variants among 12 of 71 individuals (16.9%). Among the mutation-positive subjects, 50% were diagnosed with breast cancer and had mutations in BRCA1, CDH1 and MUTYH. Notably, 33.3% were individuals diagnosed with polyposis or who had family cases and harbored pathogenic mutations in APC and MUTYH. The remaining individuals (16.7%) were gastric cancer patients with pathogenic variants in CDH1 and MSH2. Overall, 54 (76.05%) individuals presented at least one variant uncertain significance (VUS), totalizing 81 VUS. Of these, seven were predicted to have disease-causing potential. CONCLUSION: Overall, analysis of all these genes in NGS-panel allowed the identification not only of pathogenic variants related to hereditary cancer syndromes but also of some VUS that need further clinical and molecular investigations. The results obtained in this study had a significant impact on patients and their relatives since it allowed genetic counselling and personalized management decisions.

Exome Sequencing of Native Populations From the Amazon Reveals Patterns on the Peopling of South America.

Studies on the peopling of South America have been limited by the paucity of sequence data from Native Americans, especially from the east part of the Amazon region. Here, we investigate the whole exome variation from 58 Native American individuals (eight different populations) from the Amazon region and draw insights into the peopling of South America. By using the sequence data generated here together with data from the public domain, we confirmed a strong genetic distinction between Andean and Amazonian populations. By testing distinct demographic models, our analysis supports a scenario of South America occupation that involves migrations along the Pacific and Atlantic coasts. Occupation of the southeast part of South America would involve migrations from the north, rather than from the west of the continent.

Role of miRNAs in Sigmoid Colon Cancer: A Search for Potential Biomarkers.

The aberrant expression of microRNAs in known to play a crucial role in carcinogenesis. Here, we evaluated the miRNA expression profile of sigmoid colon cancer (SCC) compared to adjacent-to-tumor (ADJ) and sigmoid colon healthy (SCH) tissues obtained from colon biopsy extracted from Brazilian patients. Comparisons were performed between each group separately, considering as significant p-values < 0.05 and |Log2(Fold-Change)| > 2. We found 20 differentially expressed miRNAs (DEmiRNAs) in all comparisons, two of which were shared between SCC vs. ADJ and SCC vs. SCH. We used miRTarBase, and miRTargetLink to identify target-genes of the differentially expressed miRNAs, and DAVID and REACTOME databases for gene enrichment analysis. We also used TCGA and GTEx databases to build miRNA-gene regulatory networks and check for the reproducibility in our results. As findings, in addition to previously known miRNAs associated with colorectal cancer, we identified three potential novel biomarkers. We showed that the three types of colon tissue could be clearly distinguished using a panel composed by the 20 DEmiRNAs. Additionally, we found enriched pathways related to the carcinogenic process in which miRNA could be involved, indicating that adjacent-to-tumor tissues may be already altered and cannot be considered as healthy tissues. Overall, we expect that these findings may help in the search for biomarkers to prevent cancer progression or, at least, allow its early detection, however, more studies are needed to confirm our results.

Global Analyses of Expressed Piwi-Interacting RNAs in Gastric Cancer.

Gastric cancer (GC) represents a notable amount of morbidity and mortality worldwide. Understanding the molecular basis of CG will offer insight into its pathogenesis in an attempt to identify new molecular biomarkers to early diagnose this disease. Therefore, studies involving small non-coding RNAs have been widely explored. Among these, PIWI-interacting RNAs (piRNAs) are an emergent class that can play important roles in carcinogenesis. In this study, small-RNA sequencing was used to identify the global piRNAs expression profile (piRNome) of gastric cancer patients. We found 698 piRNAs in gastric tissues, 14 of which were differentially expressed (DE) between gastric cancer (GC), adjacent to gastric cancer (ADJ), and non-cancer tissues (NC). Moreover, three of these DE piRNAs (piR-48966*, piR-49145, piR-31335*) were differently expressed in both GC and ADJ samples in comparison to NC samples, indicating that the tumor-adjacent tissue was molecularly altered and should not be considered as a normal control. These three piRNAs are potential risk biomarkers for GC, especially piR-48966* and piR-31335*. Furthermore, an in-silico search for mRNAs targeted by the differentially expressed piRNAs revealed that these piRNAs may regulate genes that participate in cancer-related pathways, suggesting that these small non-coding RNAs may be directly and indirectly involved in gastric carcinogenesis.

Identification of NUDT15 gene variants in Amazonian Amerindians and admixed individuals from northern Brazil.

INTRODUCTION: The nudix hydrolase 15 (NUDT15) gene acts in the metabolism of thiopurine, by catabolizing its active metabolite thioguanosine triphosphate into its inactivated form, thioguanosine monophosphate. The frequency of alternative NUDT15 alleles, in particular those that cause a drastic loss of gene function, varies widely among geographically distinct populations. In the general population of northern Brazilian, high toxicity rates (65%) have been recorded in patients treated with the standard protocol for acute lymphoblastic leukemia, which involves thiopurine-based drugs. The present study characterized the molecular profile of the coding region of the NUDT15 gene in two groups, non-admixed Amerindians and admixed individuals from the Amazon region of northern Brazil. METHODS: The entire NUDT15 gene was sequenced in 64 Amerindians from 12 Amazonian groups and 82 admixed individuals from northern Brazil. The DNA was extracted using phenol-chloroform. The exome libraries were prepared using the Nextera Rapid Capture Exome (Illumina) and SureSelect Human All Exon V6 (Agilent) kits. The allelic variants were annotated in the ViVa® (Viewer of Variants) software. RESULTS: Four NUDT15 variants were identified: rs374594155, rs1272632214, rs147390019, andrs116855232. The variants rs1272632214 and rs116855232 were in complete linkage disequilibrium, and were assigned to the NUDT15*2 genotype. These variants had high frequencies in both our study populations in comparison with other populations catalogued in the 1000 Genomes database. We also identified the NUDT15*4 haplotype in our study populations, at frequencies similar to those reported in other populations from around the world. CONCLUSION: Our findings indicate that Amerindian and admixed populations from northern Brazil have high frequencies of the NUDT15 haplotypes that alter the metabolism profile of thiopurines.

NFATC2 Modulates Radiation Sensitivity in Dermal Fibroblasts From Patients With Severe Side Effects of Radiotherapy.

Although it is well established that 5 to 15% of radiotherapy patients exhibit severe side-effects in non-cancerous tissues, the molecular mechanisms involved are still poorly known, and the links between cellular and tissue radiosensitivity are still debated. We here studied fibroblasts from non-irradiated skin of patients with severe sequelae of radiotherapy, to determine whether specific basal cell activities might be involved in susceptibility to side-effects in normal tissues. Compared to control cells, patient fibroblasts exhibited higher radiosensitivity together with defects in DNA repair. Transcriptome profiling of dermal fibroblasts from 16 radiotherapy patients with severe side-effects and 8 healthy individuals identified 540 genes specifically deregulated in the patients. Nuclear factor of activated T cells 2 (NFATC2) was the most differentially expressed gene, poorly expressed at both transcript and protein level, whereas the NFATC2 gene region was hypermethylated. Furthermore, NFATC2 expression correlated with cell survival after irradiation. Finally, silencing NFATC2 in normal cells by RNA interference led to increased cellular radiosensitivity and defects in DNA repair. This study demonstrates that patients with clinical hypersensitivity also exhibit intrinsic cellular radiosensitivity in their normal skin cells. It further reveals a new role for NFATC2 as a potential regulator of cellular sensitivity to ionizing radiation.

Whole mitochondrial genome sequencing highlights mitochondrial impact in gastric cancer.

Mitochondria are organelles that perform major roles in cellular operation. Thus, alterations in mitochondrial genome (mtGenome) may lead to mitochondrial dysfunction and cellular deregulation, influencing carcinogenesis. Gastric cancer (GC) is one of the most incident and mortal types of cancer in Brazil, particularly in the Amazon region. Here, we sequenced and compared the whole mtGenome extracted from FFPE tissue samples of GC patients (tumor and internal control - IC) and cancer-free individuals (external control - EC) from this region. We found 3-fold more variants and up to 9-fold more heteroplasmic regions in tumor when compared to paired IC samples. Moreover, tumor presented more heteroplasmic variants when compared to EC, while IC and EC showed no significant difference when compared to each other. Tumor also presented substantially more variants in the following regions: MT-RNR1, MT-ND5, MT-ND4, MT-ND2, MT-DLOOP1 and MT-CO1. In addition, our haplogroup results indicate an association of Native American ancestry (particularly haplogroup C) to gastric cancer development. To the best of our knowledge, this is the first study to sequence the whole mtGenome from FFPE samples and to apply mtGenome analysis in association to GC in Brazil.

miRNome Reveals New Insights Into the Molecular Biology of Field Cancerization in Gastric Cancer.

BACKGROUND: MicroRNAs (miRNAs) play an important role in gastric carcinogenesis and have been associated with gastric field cancerization; however, their role is not fully understood in this process. We performed the miRNome sequencing of non-cancerous, adjacent to tumor and gastric cancer samples to understand the involvement of these small RNAs in gastric field cancerization. METHODS: We analyzed samples of patients without cancer as control (non-cancerous gastric samples) and adjacent to cancer and gastric cancer paired samples, and considered miRNAs with |log2(fold change)| > 2 and Padj < 0.05 to be statistically significant. The identification of target genes, functional analysis and enrichment in KEGG pathways were realized in the TargetCompare, miRTargetLink, and DAVID tools. We also performed receiver operating characteristic (ROC) curves and miRNAs that had an AUC > 0.85 were considered to be potential biomarkers. RESULTS: We found 14 miRNAs exclusively deregulated in gastric cancer, of which six have potential diagnostic value for advanced disease. Nine miRNAs with known tumor suppressor activities (TS-miRs) were deregulated exclusively in adjacent tissue. Of these, five have potential diagnostic value for the early stages of gastric cancer. Functional analysis of these TS-miRs revealed that they regulate important cellular signaling pathways (PI3K-Akt, HIF-1, Ras, Rap1, ErbB, and MAPK signaling pathways), that are involved in gastric carcinogenesis. Seven miRNAs were differentially expressed in both gastric cancer and adjacent regarding to non-cancerous tissues; among them, hsa-miR-200a-3p and hsa-miR-873-5p have potential diagnostic value for early and advanced stages of the disease. Only hsa-miR-196a-5p was differentially expressed between adjacent to cancer and gastric cancer tissues. In addition, the other miRNAs identified in this study were not differentially expressed between adjacent to cancer and gastric cancer, suggesting that these tissues are very similar and that share these molecular changes. CONCLUSION: Our results show that gastric cancer and adjacent tissues have a similar miRNA expression profile, indicating that studied miRNAs are intimately associated with field cancerization in gastric cancer. The overexpression of TS-miRs in adjacent tissues may be a barrier against tumorigenesis within these pre-cancerous conditions prior to the eventual formation or relapse of a tumor. Additionally, these miRNAs have a great accuracy in discriminating non-cancerous from adjacent to tumor and cancer tissues and can be potentially useful as biomarkers for gastric cancer.

Global miRNA expression profile reveals novel molecular players in aneurysmal subarachnoid haemorrhage.

The molecular mechanisms behind aneurysmal subarachnoid haemorrhage (aSAH) are still poorly understood. Expression patterns of miRNAs may help elucidate the post-transcriptional gene expression in aSAH. Here, we evaluate the global miRNAs expression profile (miRnome) of patients with aSAH to identify potential biomarkers. We collected 33 peripheral blood samples (27 patients with cerebral aneurysm, collected 7 to 10 days after the haemorrhage, when usually is the cerebral vasospasm risk peak, and six controls). Then, were performed small RNA sequencing using an Illumina Next Generation Sequencing (NGS) platform. Differential expression analysis identified eight differentially expressed miRNAs. Among them, three were identified being up-regulated, and five down-regulated. miR-486-5p was the most abundant expressed and is associated with poor neurological admission status. In silico miRNA gene target prediction showed 148 genes associated with at least two differentially expressed miRNAs. Among these, THBS1 and VEGFA, known to be related to thrombospondin and vascular endothelial growth factor. Moreover, MYC gene was found to be regulated by four miRNAs, suggesting an important role in aneurysmal subarachnoid haemorrhage. Additionally, 15 novel miRNAs were predicted being expressed only in aSAH, suggesting possible involvement in aneurysm pathogenesis. These findings may help the identification of novel biomarkers of clinical interest.

The comprehensive expression analysis of circular RNAs in gastric cancer and its association with field cancerization.

Circular RNAs comprise a new class of long noncoding RNAs characterized by their 5' and 3' ends covalently joined. Previous studies have demonstrated that some circular RNAs act as microRNA sponges, and are associated with cellular proliferation in cancer. We were the first to analyze the global expression of circular RNAs in samples of patients without gastric cancer, gastric cancer, and matched tumor-adjacent gastric tissue. Among the samples, we identified 736 previously annotated circular RNAs by RNA-Seq. The tumor-adjacent tissue presented the higher abundance of circular RNAs and could not be considered as a normal tissue, reinforcing the notion of field effect in gastric cancer. We identified five differentially expressed circular RNAs that may be potential biomarkers of this type of cancer. We also predicted candidate microRNAs targets of the highest expressed circular RNAs in gastric tissues and found five miRNAs. Overall, our results support the hypothesis of circular RNAs representing a novel factor in the dynamic epigenetic network of gene regulation, which involves the microRNAs, its mRNAs targets, and the circular RNAs-derived genes. Further studies are needed to elucidate the roles and the functional relevance of the circular RNAs in human diseases.